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| We try to understand flavor oscillations and to develop the formulae for describing 

& ■ neutrino oscillations in de Sitter space-time. First, the covariant Dirac equation is 

investigated under the conformally flat coordinates of de Sitter geometry. Then, we 
obtain the exact solutions of the Dirac equation and indicate the explicit form of 
the phase of wave function. Next, the concise formulae for calculating the neutrino 
oscillation probabilities in de Sitter space-time are given. Finally, The difference 
I/-) ■ between our formulae and the standard result in Minkowski space-time is pointed 

i out. 
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1 Introduction 



Neutrino oscillations have been the most natural explanation to the so-called 
solar neutrino problem and atmospheric neutrino problem [1]. Since neutrino 
oscillations may provide a tool both to understand elementary particles and 
their interactions and to uncover the secrets on ancient celestial bodies, they 
have attracted a lot of attention in modern physics [2]. In Minkowski space- 
time, neutrino oscillations have been extensively studied in the past by using 
both plane waves [3] and wave packets [4,5] to represent the emitted neutrinos. 
But, the neutrinos emitted by ancient sources must travel a very long distance 
to arrive at our detectors, and be affected by gravitation. Therefore, in more 
recent years, physicists have turned their attention to specifically gravitational 
contributions to neutrino oscillations. 
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According to the methods appeared in the literature, we divide those con- 
siderations into three categories: 

(1) Using approximate operator [6,7,8,9]. In these papers, the effective Hamil- 
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tonian or the semiclassical approximation, such as, p x ~ E — is gen- 
erally used to investigate the complicated physical system. This kind of 
method is fit for considering matter effects on neutrino oscillations. 

(2) Adopting semiclassical phase factor [10,11,12,13]. This category of study 
first appeared in Ref. [10]. All of these papers laid their foundation on 
an extended covariant phase factor, for instance, $ = p^dx^, given in 
Ref. [14]. Since the semiclassical action is taken as the quantum phase, 
couplings between spin and gravitation are thus ignored. 

(3) Studying the covariant Dirac equation [15]. It seems that this is the unique 
way to acquire the exact formulae for calculating the oscillation probabil- 
ities. But, only in several ideal models, one can obtain the exact solutions 
of the Dirac equation. So this method can be used in a few cases. 

Recent astronomical observations on Type la supernova [16,17,18] and the 
cosmological microwave background radiation [19,20,21] imply that dark en- 
ergy can not be zero, and which will dominate the doom of our Universe. 
Therefore, our observable universe will be an asymptotical de Sitter space- 
time. On the other hand, the scenarios of inflation usually need an inflation- 
ary epoch, in which space-time can be treated as a de Sitter manifold [22]. 
Upon above considerations, we investigate neutrino oscillations in de Sitter 
space-time in this paper. 

The different evolution of the phases of mass eigenstates plays the pivotal 
role in neutrino oscillations. In Minkowski space-time, the phase of a free boson 
or of a free fermion is given by solving the Klein-Gordon equation or Dirac 
equation, that is, x) = —{E -t — p-x). From this point of view, we need 
to solve the covariant Dirac equation to get the exact form of the phase of 
wave function. After that, we can obtain the exact formulae for calculating 
the neutrino oscillation probabilities. 

This paper is organized as follows. The covariant Dirac equation in de Sitter 
manifold is explicitly given in section 2. In section 3, the exact solutions of 
the Dirac equation are obtained. After that, we read out the phase factor. 
We then present the formulae on neutrino-oscillation probabilities in section 
4. Conclusions and remarks appear in the last section. 
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2 Dirac equation in de Sitter space-time 



To set the stage for studying neutrino oscillations in a curved space-time, 
we introduce the local tetrad field as follows 1 

gT = e^(x)e\(x) V ab , (1) 



where g^ u being the space-time metric, and the metric tensor of Minkowski 
space-time r\ ah satisfy 



rf = +1 , r] 11 =71™ =7) 



11 _ ^,22 _ 33 



r] ab = for a ^ b 



(2) 



In a curved space-time, the covariant Dirac equation reads 2 

" d 



dx^ 



r M (x) 



ip(x) = m'4>(x) 



(3) 



Here 7 M (x) are the curvature-dependent Dirac matrices and r M (x) are the spin 
connections to be determined. The curvature-dependent Dirac matrices ^{x) 
are presented in terms of the tetrad field as 



-f(x) = e» a {x)r 



(4) 



where 7° denotes the standard flat-space Dirac matrices, which satisfies {7", •y b } = 
2r) ab . 

Under a general coordinate transformation, x M — >• y u , it is well known that 
the spinor keeps invariant [23], that is, 

^(y) = ^0) • ( 5 ) 



Hence the phase difference of two points for a free neutral spin-| particle is 
invariant under the coordinate transformation. As we will demonstrate, the 
phenomenon of flavor oscillation is only related with the phase difference and 
the mixing angles, one then can select any convenient coordinates to study 
flavor oscillation. Because of this property, we choose the comoving coordinates 
to investigate the neutrino oscillations in de Sitter space-time. 

1 We use Roman suffixes to refer to the bases of local Minkowski frame, and use 
Greek suffixes to refer to curvilinear coordinates of space-time. 

2 we adopt the conventions that x° = t, x 1 = r, x 2 = 9, x 3 = cp and h = c = 1. 
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In the conformally flat coordinates the metric of de Sitter space-time is 
taken to be 



ds 2 = g^dxTdx" = dt 2 - R 2 (t) \dr 2 + r 2 (d9 2 + sin 2 



(6) 



where R(t) = e Ht is called cosmic scale factor of de Sitter space-time, if is a 
constant. 

The tetrad field in the case of de Sitter metric (6) is thus taken to be 



e>M a( x ) = dia S ( 1> 4> 7T' 



R Rr Rr sin 9 



(7) 



Inserting Eq.(7) into Eq.(4) yields 



7°(a;)=7o 



7 (a?) = -^7i 



Rr sin # 



^^73 • 



The spin connections T At (x) in Eq.(3) satisfy the equation 



(9) 



where are the Christoffel symbols for the metric (6), which are determined 
by 



r = -s VT ( ^ + ^ - dgw 

w 2 g I dx* dxP dx T 



(10) 



Hence Eq.(9) can be solved for the spin connections T^x) which we determine 

as 



r = o 



r R 

1 1 = -^iQll , 



1 ■ 1 

T 2 = --Rr7o7 2 + -727i , 

1 • 1 1 

T 3 = -Rr sin ^7073 + - sin 6*737! + - cos 6*7372 
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where R denotes so that the combination 7 At (x)r /J (x) in Eq.(3) simplifies 
to 



r{x)VJx) 



3R 



cos 9 



2R l0+ Rr 11 + 2Rrsin9 



72 



(12) 



Inserting Eq.(8) and Eq.(12) into Eq.(3), then we obtain the Dirac equation 
in de Sitter space-time 



1 d 
H ^t73- 



d l\ 1 ( d cos 6 



7ihr + - +-72U 



r sin 9 dip 



dr r I r " \89 2 sin 9 



(13) 



3 The exact solutions of the Dirac equation 

We can simplify the Dirac equation (13) by transforming [24,25] 

i/j(t,r,9,ip) = (fBn6)-sR-*V(t,r,9,<p) . (14) 



The reduced equation in terms of ty(t, r, 9, ip) is of the form 



7071 H + I) + J 7li ^'^ 



(15) 



where 



j ^ )=7 ° 71 ( 7 4 +73 ^?^ 



(16) 



is a Hermitian operator, as we know, which is defined in Ref. [26] first. The 
eigen equation of K is 



KV^ = q% tK , ? = 0,±1,±2, 



(17) 



where the eigen value q is an integer, which has been proven in Ref. [23,26]. 
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Only the two matrices 70 and 71 remain explicitly in the simplified Dirac 
equation (15), they can therefore be represented by 2 x 2 matrices 



7o 



1 \ / -n 

0-1 In 



where 1, n are respectively 2x2 unit matrix and the Pauli matrix, namely 

ri = I ° 1 I • (19) 
1 



We separate the angular factor from the wave function, and represent the 
radial and temporal factor by a two-component spinor, that is 

*„ = e<{o,tp) I 0M) I , (20) 



where the angular factor, ; (#, <p), is determined by the requirement KQ q (6, ip) 
t;O f (0, ip). This equation for eigenstates of the angular motion has been inves- 
tigated by Schrodinger [26]. He found, the operator K is related to the total 
angular momentum. Their eigenvalues satisfy [27] 

l f-(/ + l) for 3 = 1 + 1- 
? = =F(j + 5 )= \ ■ (21) 

[l for j = I - - 



Also the eigenfunction f is related with the spherical harmonics as follows 
B q (9, ip) oc Y m+ 1 (6, <p) oc Pr +I (cos 6)e^ + ^ , (22) 

or 

9,(0, <p) oc Y hw _i{0,<p) oc Pr^(cose)e^^^ , (23) 
where P(cos 6) being the associated Legendre functions, and m±| = 0, ±1, ±2, 

;±l. 
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With the help of equations (17), (18), and (20), we can rewrite the equation 
(15) into two equations 



R + irnj <f>(r,t) 



d_ r 

dr r 

'd_ I s 

dr r 



nx(r,t) , 

Tx(f)(r,t) . 



(24) 
(25) 



Here -T\q is the term of spin-orbit coupling. To solve the above equations, 
we separate the functions 0(t, r) and x(t, r) into radial and temporal factors, 
respectively 



<f>(t, r) = U^T^t) , X (t, r) = U 2 {r)T 2 {t) 



(26) 



With the help of Eq.(19), inserting the above equations into Eq.(24) and 
Eq.(25) yields the evolution equations 



R ( — + im J Ti(t) - iKT 2 (t) = 



R\ — -im \ T 2 {t) -wT x {t) = 



and the radial equations 



(27) 
(28) 



' d_ r 

^dr r / 

' &_ 1 N 

dr r , 



' d_ r 

dr r 



dr r I r 



U x + k U i = 
U 2 + k 2 U 2 = 



(29) 
(30) 



Obviously, if U\ is substituted by U 2 , and simultaneously <^ are replaced by —q, 
then the equation (29) becomes the equation (30). Therefore, in the following, 
it is of only necessity to study the equation (29). 



The equation (29) directly simplifies to [25] 



d 2 U x 2 dUt 



dr 2 



+ 



r dr 



k ) U = 



(31) 



Through a simple transformation, one can easily find that the equation (31) 
is the Bessel's equation of order —q + ^. The solutions of the above equation 
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are [28] 



UAr) 



-.J. 



'nr 



Kr 



(32) 



Where J±^-i)(Kr) are called the Bessel functions. In the case of q = 0, —1, —2, ■■ 
•, the Bessel functions in terms of sine or cosine functions are explicitly taken 
to be [28] 



^r\ l ^ J (_i)"(_ ? + 2n)! 



2 V 7 V mvr I V 2 / ^ (2n)!(-^-2n)!(2Kr) 2n 



(33) 

+ cos I Kr + — J ^ 

Z 7 71=0 



.s±li 

2 



(2n + !)!(-? - 2n - l)!(2Kr) 2 ™+ 1 



and 



J ? _i (Kr) = V — S cos I «r - — 2^ 



2 V 7 V 7T«r j V 2 / ^ (2n) !(-?- 2n)!(2Kr) 2 ™ 

(34) 

(-l) n (-q + 2n + 1)! 



sm Kr 



V 

2 ^ ^ ( 2ra + 1)K-? -2n- l)!(2Kr) 2n+1 



The formulae for <j > are easy to read out from above two equations. 

We have indicated that the solution of can be obtained from the solution 
(32) by substituting — q for q. Hence we have acquired the exact forms of the 
spatial factor of wave function in de Sitter geometry. 



4 The factor of time evolution 



After setting s = 1/R(t), since R(t) = e Ht , it is easy to get the following 
formula 

ff=i^ = -!£. (35) 

R dt sdt v ' 



In terms of new variable s, the evolution equations (27) and (28) can be 
rewritten as 
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im. 



- HT[ + — 7\ - ikT 2 = , 
s 



im 



ET' 2 T 2 - mTi = , 



(36) 
(37) 



where T' = dT/ds. From the above equations, we acquire the equation satisfied 
by T\ as follows 



s 2 T'l + 



H 2 



v{v-\) 



Ti = 0, 



(3* 



where 



im im 

_ or r/ = 1 - — . (39) 



Obviously, the above equation is the Bessel's equation of order v — |, the 
solutions of which is generally presented as 



These solutions demonstrate that v — ^ and z/ = 1— ^ are quite equivalent. If 
2u isn't an integer, then the Bessel function J±^_i^ can be explicitly expressed 
as the infinite series, that is 



j ( JL.) = V 1 f — Y +U ~ l 

V -^HR> ^ t! T(v + l + \) \2HRJ 

where T(- • •) is the gamma function. It is easy to read out the form of T 2 (t) 
T 2 (t) = -^=J +( ,m m(— ) • (42) 



Hence we have given the exact solutions of the Dirac equation in de Sitter 
space-time. 

In a different coordinates of spatially flat Robertson- Walker space-time, the 
exact solutions of the Dirac equation with R(t) = e Ht have been given in 
Ref.[29]. We acquire the same evolution factor as the result in Ref.[29]. 
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5 Neutrino oscillations 



Recent neutrino experiments imply that the masses of neutrinos satisfy 
m ~ 10~ 2 eV [2]. If one select H as the present-day Hubble constant, then 
\v\ is about 10 30 . Therefore, \v\ is a huge number. Then we can acquire the 
approximate formulae of the Bessel function from Eq.(41) 




where c\ and C2 are small constants. From above equation, we can reasonably 
obtain more concise expression as follows 



K 



HR' 



exp < v 



v 



log 



( K 



\2im 



I im 
1 ex P i -Jj 



log 



Ht 



(-) 

\2m) 



imt 



(44) 



It is well known that the Dirac equation describes both particle and its an- 
tiparticle. So the sine or cosine functions shall be rewritten as the exponential 
function with imaginary argument. According to the exact solutions of the 
Dirac equation, we can read out the general form of the phase in de Sitter 
space-time, that is 



m 

phase = — 
H 



1 + log — 

2m 



1 ^7T 

mt + nr + (m ± -)ip ± — + z^O) , (45) 



where z^iff) is a function related with eigenvalues q, m and coordinate 9. 

Updated large experiments have shown that neutrinos are possibly massive 
although their masses are very tiny. If there are flavor mixing in neutrino 
sector, and the leptonic quantum numbers are only approximately conserved, 
then there exist the fascinating possibility of neutrino oscillations. 

In the case of three generations of neutrinos, each of the three neutrino mass 
eigenstates shall be represented by Vj\ I = 1,2, 3. Because of neutrino mixing, 
we have the linear superposition 

vk= E tfAifi, (46) 

7=1,2,3 
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where Uai is a unitary matrix, and A = e, [i,t represents the weak flavor 
eigenstates (corresponding to electron, muon and tau neutrinos respectively). 
The Ui, v 2 and u 3 correspond to the three mass eigenstates of masses mi, m 2 
and 7713, respectively. Conventionally the mixing matrix U A i is parameterized 
by three mixing angles, one CP violating phase and two Majorana phases. For 
the sake of investigation on neutrino oscillations, we ignore the CP violating 
phase. In this paper we assume that the neutrinos are of the Dirac type. Hence, 
in terms of mixing angles a, /3, rj, the mixing matrix reads [2] 



U(a,P,rj) 



-C a S/3S v - S a Cf3 +C a Cp - S a SpS v SpC^ 
y -CaCpSr, + S a S{3 -C a S{3 - S a CpS rj CpC^ J 



(47) 



where = cos(£), = sin(^), with £ = a,[3,r). 

Let us assume: a weak flavor eigenstate, e.g., an electron neurtino, is emitted 
at a point P with coordinates P(t e , r e , 9 e , <p e ). For simplicity, we are interested 
in its radial propagation. We also assume that the different mass eigenstates 
have the same total angular momentum So, in the momentum representa- 
tion, the mass eigenstates must have different eigenvalue k, denoted by kj. 
Therefore, in de Sitter space-time, the radial evolution of them is given by the 
expression: 



uj(t > t e , r) = exp 



f .777/ 


1 + log 




Vir 




V2m,/J 



-imi(t - t e ) + i«/(r - r e ) + iU} uj(P) , (48) 

here U being the phase factor which is the same for all mass eigenstates. 

Considering Eq.(47) and Eq.(48) together, one can easily draw the following 
conclusion: in de Sitter space-time, under the assumption that the mass eigen- 
states have the same total angular momentum, if a weak flavor eigenstate is 
created at the point P, then the probability to find the weak flavor eigenstate 
v\i at the point Q(tf, r/, 9 e , y? e ) is given by the formulae: 



p(u A -> v A ,) =| u A ,(Q) | 2 

3 

= {UA>iexp[-imi(tf -t e ) +iK I (r f - r e )]U AI } 
i,j=i 

■ {U A/J exp [+imj(t f - t e ) - inj(r f - r e )\ U AJ ) 
= 5 A , A - ^U A 'iU A1 U A > 2 U A2 sin 2 [u;i2] - 4f/ A 'i^Ai^A'3^A3 sin 2 [^i 3 ] 
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- ^U a/2 Ua2U a/3 Ua 3 sin 2 [cu 23 ] , (49) 

where 



Ar 
Ak 



Ak a 

Ar 

2 



Am 



At 



Am 



log 

2H 6 



' e , 



2# 

At = t f -t e , 
Am = mj — mi 



f kj 
[\2mj 



mj 



(2m\ 



\ Ki 



m i 



(50) 

(51) 
(52) 



We have given the formulae for calculating the neutrino oscillation proba- 
bilities in de Sitter space-time, which demonstrate that the probabilities are 
related with not only the positions, masses and momenta, but also the Hubble 
constant H. In Minkowski space-time, the Hubble constant doesn't appear. 



6 Conclusion 



The neutrino oscillations in de Sitter space-time is studied through investi- 
gating the covariant Dirac equation. The exact solutions of the Dirac equation 
are obtained in the conformally fiat coordinates. The phase of the spin-^ wave 
function is read out. Then the formulae of neutrino oscillation probabilities 
in de Sitter space-time are constructed. Our result is quite different from the 
stand form in Minkowski space-time because of the different topology of de 
Sitter manifold. 
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